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Abstract 

The theoretical description of strongly correlated quantum systems out of equilibrium 
presents several challenges and a number of open questions persist. In this paper we focus 
on nonlinear electronic transport through an interacting quantum dot maintained at finite 
bias using a concept introduced by Hershfield [Phys. Rev. Lett. 70, 2134 (1993)] whereby 
one can express such nonequilibrium quantum impurity models in terms of the system's 
Lippmann-Schwinger operators. These scattering operators allow one to reformulate the 
nonequilibrium problem as an effective equilibrium problem associated with a modified 
Hamiltonian. In this paper we provide a pedagogical analysis of the core concepts of the 
effective equilibrium theory. First, we demonstrate the equivalence between observables 
computed using the Schwinger-Keldysh framework and the effective equilibrium approach, 
and relate the Green's functions in the two theoretical frameworks. Second, we expound 
some applications of this method in the context of interacting quantum impurity models. We 
introduce a novel framework to treat effects of interactions perturbatively while capturing 
the entire dependence on the bias voltage. For the sake of concreteness, we employ the 
Anderson model as a prototype for this scheme. Working at the paiticle-hole symmetric 
point, we investigate the fate of the Abrikosov-Suhl resonance as a function of bias voltage 
and magnetic field. 
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1 Introduction 



Quantum systems exhibiting an interplay of interactions and out-of-equilibrium ef- 
fects are of significant interest and constitute an active area of research | 
In contrast to the equilibrium setting, the theoretical foundations of these systems 
are yet to be well understood. Questions such as whether renormalization group 
and universality concepts can be generalized to the nonequilibrium domain persist. 
Furthermore, there is a lack of a unifying framework for treating such systems. 
The real-time Schwinger-Keldysh formalism II11I12I13I14II has been successfully 
employed to study certain classes of problems. Most such treatments have been 
perturbative and is often a delicate procedure. Care must be taken to avoid the ap- 
pearance of infrared divergences, an artifact of the the infinite contour of integration 
[[T5l . While quantum impurity models in equilibrium have been extensively studied 
[[T6l , probing and understanding properties of nonequilibrium steady states is a far 
more subtle task. Nevertheless, there have been significant advances in our under- 
standing via several distinct approaches, which include the scattering Bethe Ansatz 
[I17I18II , field theory techniques II19I20II , time-dependent density matrix renormal- 
ization group (RG) [i20i21J . time-dependent numerical RG II22I23II . perturbative 
RG 1124125 126H . Hamiltonian flow equations (23, functional RG 112812913011 . strong- 
coupling expansions B1I32I33I34II . master equations (351, diagrammatic Monte 
Carlo B6I37I38I39I40II and imaginary-time nonequilibrium quantum Monte Carlo 

mm. 



Mesoscopic quantum systems such as quantum dots ("artificial atoms"), which are 
characterized by a set of discrete energy levels, when coupled to reservoir leads (43l 
are described by quantum impurity models. In this context, the Anderson model, 
which describes a single discrete level coupled via tunneling to a Fermi-liquid sea, 
is of particular interest (T6l . For a level with sufficiently low energy and with strong 
on-site interaction, the Anderson model mimics a Coulomb-blockaded quantum dot 
(441 . In this regime, charge degrees of freedom on the level are frozen out and the 
Anderson model becomes intimately related to the Kondo model. The latter de- 
scribes a magnetic impurity (manifested by the spin of the highest occupied level 
for an odd number of electrons on the quantum dot) entangled with the spins of 
a Fermi sea (In the mesoscopic setting, the reservoir leads play the role of Fermi 
seas). The Kondo entanglement then produces a prominent Abrikosov-Suhl reso- 
nance in the density of states of the quantum dot at the Fermi level, which results in 
perfect transparency when the quantum dot is coupled symmetrically to the source 
and the drain II43I45I46I47I48L Nanoscale systems can be routinely driven out of 
equilibrium by applying a bias voltage between the two reservoir leads. Among 
the various issues which arise out of equilibrium, the precise fate of the Abrikosov- 
Suhl resonance (49l when applying a finite bias voltage, remains an issue of debate. 
In general, it seems essential to elaborate theoretical and numerical methods which 
will allow access to the full current-voltage characteristics of interacting meso- 
scopic systems. 
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The main purpose of this paper is to reformulate electronic nonequilibrium trans- 
port in quantum impurity models in terms of an effective equilibrium steady- state 
density operator. This concept was initiated by Hershfield ||50l who proposed that 
the appropriate density operator can be expressed in terms of the system's Lippmann- 
Schwinger operators in a Boltzmann form. While noninteracting and effectively 
noninteracting theories lISTI have been treated analytically within this framework, 
extensions to interacting models have yet to be explicitly demonstrated. Recently, 
there have been several numerical efforts at implementing this scheme II22I23I41I42II . 
Also, the interacting resonant level model has been treated perturbatively using the 
steady-state density matrix along with dynamical "impurity conditions" ||52l , ef- 
fectively circumventing explicit construction of Lippmann Schwinger states. Given 
the growing interest in this alternative description, it is important that the theoretical 
foundations of the effective equilibrium theory be more concretely established and 
details of the formulation refined. In this paper we provide a pedagogical exposition 
of the core concepts of the effective equilibrium theory. Furthermore, we introduce 
a novel perturbative framework for interacting systems, which makes explicit use 
of its Lippmann-Schwinger scatting states. 

First, we discuss the effective equilibrium theory by making use of the open system 
limit ||53l which guarantees that the system relaxes to a steady state. This well- 
controlled mathematical construct mimics the role of relaxation mechanisms, and 
does not require the explicit inclusion of the bath degrees of freedom at the level of 
the Hamiltonian. We also demonstrate the equivalence between observables com- 
puted using the Schwinger-Keldysh framework and the effective equilibrium ap- 
proach. It is important to emphasize the fact that the effective equilibrium descrip- 
tion encompasses the case when one includes interactions on the dot. 

Next, we outline a systematic scheme for computing Green's functions perturba- 
tively in the interaction strength, which also provides the basis for an extension to 
non-perturbative analytical methods. Having established the basic methodology, we 
use the Anderson impurity model [[54ll out of equilibrium as a prototype to illustrate 
the details of our perturbative approach. Working at the particle-hole symmetric 
point, we derive the self-energy of an electron on the dot to second order in the in- 
teraction strength, capturing the complete (non-linear) bias dependence. This result 
is then used to study the spectral function within the Born Approximation. The ob- 
served central peak can be regarded as a precursor of the much sharper Abrikosov- 
Suhl (Kondo) peak. Our results correctly reduce to the well-known Fermi liquid 
expansion in the limit of small bias, temperature and frequency [I16I55II . When in- 
creasing the bias voltage, there is a complete suppression of the Kondo effect and 
no resonance splitting is observed. We then break particle-hole symmetry by in- 
troducing a local magnetic field on the dot and analyze the effect of the bias and 
the Zeeman energy on the Abrikosov-Suhl resonance ||39l . This approximation is 
known to qualitatively reflect all features one expects in the low energy sector (in 
particular at low bias voltages we obtain a typical trident-shaped curve). 
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Our paper is organized as follows. In Section 2, we introduce the model and define 
the "open-system limit" I118I53B . This concept provides a transparent framework 
that guarantees the existence and uniqueness of a steady-state density matrix. In 
Section 3, we summarize the effective equilibrium approach introduced by Hersh- 
field [50J and demonstrate the equivalence between observables computed within 
this framework and the Schwinger-Keldysh formalism. This enables us to arrive 
at the Meir-Wingreen formula for the current I156I57II . We then derive a compact 
expression for the charge occupation on the dot, whereby we identify the spectral 
function of the dot as the central quantity of interest. In Section 4, we introduce 
a novel perturbative framework for treating interactions which captures the effect 
of the bias voltage non-perturbatively. In Section 5, we present our results for the 
current and charge occupancy on the dot for the Anderson impurity model ll54ll . 
The details of the calculations are included in the appendices. 



2 Model 



In our discussion of the effective equilibrium approach, we consider, for the sake 
of concreteness, a system consisting of two Fermi-liquid leads, coupled by tunnel 
junctions to a central system with a single energy level, the "dot" (see Fig. [T]). The 
Hamiltonian of such a generic system can be written in the form H = + + 
Ht- The first term 

Hl = Y1 (^akc\ka(^aka (1) 
akcr 

describes the left and right leads {a = ±1), where Cako- (.cj^ka) annihilates (creates) 
an electron (strictly speaking a Fermi-liquid quasiparticle) in state k with spin pro- 
jection cr in lead a. The corresponding energy dispersion is denoted by eak- These 
leads couple to the dot, whose Hamiltonian is given by 

Hd= Y.^ddid^ + ^inu (2) 

cr 

where (cij.) annihilates (creates) an electron with spin a in the discrete level with 
energy e^. Any electron interactions are lumped into the contribution Hint, and we 
will assume that these interactions are localized on the dot (as it is the case for 
the Anderson model). Finally, the tunneling of electrons between leads and dot is 
captured by the tunneling Hamiltonian 

HT = ^T.tak{cikJ. + h.c.), (3) 

V" aka 

where is the lead volume (assumed identical for both leads) and tak specifies the 
tunneling matrix element for electron transfer between state k in lead a and the 
discrete dot state. In the presence of a bias voltage, realized as a chemical potential 
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difference $ = /ii — between left and right lead, the tunneling induces an 
electric current. 

Typically, the steady state is reached after a short time determined by relevant re- 
laxation rates in the leads. For the purpose of calculations aiming at steady-state 
quantities, it is convenient to avoid the consideration of microscopic relaxation 
mechanisms and instead invoke the so-called "open-system limit" II18I53II . In short, 
this approach proceeds as follows: Initially, up to some time t = to < in the early 
past, the tunneling term is absent and the Hamiltonian of the system is given by 
Hq = Hl + Hf), such that leads and dot are decoupled. For t < to the system is 
hence described by the separable density matrix 



Here, (3 = (kBT)~^ denotes the inverse temperature and Na = Y^ka ^aka^aka is the 
number operator of electrons in lead a. Between times to < ^ < 0, the tunneling 
is then 'switched on' adiabatically, i.e., H = Hq + HTe^^9{t — to), where the 
parameter — )• 0+ defines a slow switch-on rate. At time t = the tunneling has 
reached its full strength, the system is in its steady state and observables can be 
evaluated. 

As demonstrated in Refs. IfTSl and [|53l , the existence and uniqueness of the steady 
state is tied to the validity of the inequalities vp/L <^ l^ol ^ ^ where vp de- 
notes the Fermi velocity and L the linear system size. Intuitively, these inequalities 
ensure that hot electrons hopping onto a given lead at time to will not be reflected 
back and return to the junction before the measurement process, and further that the 
process of switching on Ht remains adiabatic. We also note that the energy scale of 
switch-on |to|^^ suffices to smear out the energy level spacing vp/L. In this sense, 
the openness of the system provides the "dissipation" mechanism necessary for 
the steady state, allowing the high-energy electrons to escape to infinity and thus, 
effectively relax. On a more technical level, the inequalities result in the factoriza- 
tion of long-time correlation functions, which facilitates the proof of existence and 
uniqueness of steady state (see Appendix A and Appendix B). 

The crucial ingredient in our discussion is that L/vp determines the largest time 
scale in our problem. The exact protocol by which we switch on the tunneling is 
irrelevant for the formation of steady state. To illustrate this fact, take for example 
the non-interacting resonant level, and for convenience assume that the leads are 
identical (i.e., tak = t). Now, in stark contrast to the open-system limit we switch 
the tunneling on instantaneously. One then observes that the transients decay with 
a relaxation time ~ l/F, where F = 27rt^i/ denotes the linewidth of the dot ll58l . 
Here symbolizes the density of states, which is assumed to be constant. To make 
the argument concrete, assume that the tunneling Hamiltonian is absent for t < 
0, and the distribution functions in the leads are given by Fermi functions /(e — 




(4) 
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a*). Furthermore, suppose that the dot is initially unoccupied, i.e., = 0- The 



occupation of the dot at a time t(> 0) is given by p8 



Udit) = rid,^ (l + e-2r*) 



2Te 



-rt 



2' 



2' 



cos [(g; - €d)t] 



(5) 



Here 



du 



/(u;-f) + /(a; + f) 



(6) 



denotes the steady state expectation value of the dot occupation, after the transients 
have decayed. It is interesting to note that if we turn the coupling of the QD to the 
left and right leads to zero at the same rate, (formally this implies taking F ^> in 
Eq. ^) then the final QD occupation is given by 



^d 



fi^d-^ 



f(^d+- 



(7) 



The fact that n'^ 



n 



fin 



d ' 



clearly illustrates the irreversibility of the turning on 



process. Thus, the tunneling is treated exactly within the effective equilibrium ap- 
proach. Interactions on the other hand can be turned on adiabatically and this pro- 
cedure is completely reversible, and does not exhibit similar anomalies. It is thus 
possible to treat interactions perturbatively. 



3 Effective Equilibrium Theory 



The description of a nonequilibrium problem in terms of an effective equilibrium 
density matrix was first proposed by Hershfield in Ref. [[50| . This Section is orga- 
nized as follows. First, we outline the theoretical fundamentals of this equilibrium 
description. Then, we show the equivalence between observables computed within 
the Schwinger-Keldysh scheme and the effective equilibrium theory. In Appendices 
A and B, we flesh out the proof of the existence and uniqueness of the effective 
equilibrium density matrix, making careful use of the "open-system limit" II18I53II . 

3.1 Effective Equilibrium Density Matrix 

The nonequilibrium steady-state density matrix p is rewritten in the usual Boltz- 
mann form, 

p = exp [-^{H - Y)] , (8) 
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at the cost of introducing an additional operator, which we will call (following 
Hershfield's convention) the Y operator^ Formally, the definition of this correction 
operator as 

Y = hnp + H (9) 

is always possible. However, in this form it is neither particularly elucidating nor 
useful for calculating nonequilibrium transport properties. Hershfield put forward 
the idea that Y can be expressed explicitly and compactly in terms of Lippmann- 
Schwinger operators II59I60II i/Jaka, which are fermionic operators that diagonalize 
the full Hamiltonian [|6T1l 



H = J2eakVa^a'lpaka, (10) 
ak 

and can be formally expressed by the equation 



i'ika — r I • ^TC^aka- (H) 

eak - C + IT] 

The Liouvillian superoperators C, Ct and Cy are defined such that that their action 
on any operator O is given by CO = [H, O], CtO = [Ht, O] and CyO = [F, O] 
respectively. 

Hershfield showed that the Y operator has the general form [|50ll 

$ , 

^ = ^Z]«^afc>afc<T- (12) 

aka 

Since this operator encodes the entire $ dependence, the Y operator is also called 
the bias operator. Note that Y vanishes at zero bias and the steady-state density 
matrix correctly simplifies to the equilibrium density matrix. 



In Appendices A and B, we provide a detailed derivation of Y given in Eq. (112]) . 
In contrast to Hershfield, we do not invoke the presence of additional relaxation 
mechanisms to reach steady state. Instead, we make systematic use of the time- 
dependent open system approach [1181531 . 



3.2 Correspondence with the Schwinger-Keldysh Approach 



We now illustrate how to compute transport observables such as the current and 
the charge occupation on the dot within this description. We propose an imaginary- 
time formulation for treating such nonequilibrium systems in a manner similar to 
finite temperature field theory. 

^ Note that throughout the text, density matrices are not assumed to be normalized. Nor- 
malization is always introduced explicitly when evaluating expectation values. 
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In imaginary-time, we define the propagation of an operator by 

0{t) = e-(^-^)Oe-^(^-^) = e^(^-^^-)0. (13) 
The nonequilibrium thermal Green's function is defined on < r < /3 as 

Qo.oAr) = -(T [Oi(r)O2(0)]) = -(Oi(r)O2(0)). (14) 
The subsequent Fourier transform in imaginary-time results in 

Qo^oA^^n) = {{Oi, -. 7r—^02 ) , (15) 

where u;„ = (2n + l)vr//3 {n E Z) denotes the fermionic Matsubara frequencies. 

Switching to real-time, the Heisenberg representation of an operator O is given by 
0{t) = e*^*Oe~*^*. The nonequilibrium real-time retarded Green's function can 
then be expressed as 



tr 



By using the spectral representation and then Fourier transforming, we obtain 

We emphasize that despite the effective equilibrium character of the Hershfield ap- 
proach, one important difference between the effective equilibrium description and 
the traditional equilibrium many-body formalism remains. This distinction arises 
from the different propagators in imaginary versus real-time; namely, in Hersh- 
field's effective equilibrium formalism all Heisenberg operators in imaginary-time 
evolve under the modified Hamiltonian H — Y, whereas Heisenberg operators in 
real-time evolve under the Hamiltonian iJ. As a result, imaginary-time and real- 
time Green's function, Goi02{^^n) and G^^\q^{uj) axe not related via the analytic 
continuation iun cu + ir]. 

As an observable of prime interest in transport, let us now consider the current 
flowing through the dot. It is given by 

^ ^ h + 1^1 ^ _e / d(iVi(t)-iV_i(t)) \ 
2 2 \ dt I 

• ^ /Fat TTl\ • 6^afc 



^ a aka ^V" 

Im a—^Q ,t (t = 0) 

.aka V " 



(18) 
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where — e denotes the charge of the electron. Here, we have used the fact that the 
system is in a steady state and the current is unchanged under time translations. 
For a diagrammatic analysis, it is convenient to express the current in terms of the 
Fourier representation of the imaginary-time Green's function, namely 



Im 



E 

OckcFUJrL. 



etak 1 ^ 



t [lUJn 



(19) 



We now show the equivalence of this approach with the Schwinger-Keldysh for- 
malism and recover the familiar Meir-Wingreen formula for the steady-state current 
II56I57L For simplicity, let us assume that the coupling to the dot is independent of 
k, i.e. tak = ta, and that the leads are identical eak = ^k (this is not a strict require- 
ment and the proof can be easily generalized). Our starting point is Eq. (fT9l ) and we 
now write explicitly 

Using Eq. ([TT]) for Caka and evaluating the effect of Ct, we obtain 
1 1 




^ -d^, . . (21) 



ek + L — 17] luj. 



Making an additional observation, one can show that the second term does not 
need to be evaluated when computing the current /. To see this, we exploit the 
fact that in steady state, the current in the left and right junctions must be iden- 
tical. Hence, we can write the current as a weighted average of the form / = 
(tll-i + f^ih^ I {t\ + i^i)- This eliminates the presence of the second term in 
Eq. (I2TI) from the expression for the current. Proceeding with the remaining term 
we obtain in several steps. 



^" fU-a%\G-\,{e^), (22) 



^ \ 2/ 



where the transition from the first to the second line is facilitated by the general 
relations summarized in Appendix C. 

Using this identity together with the replacement ^ X^fc ~^ / ^ c^^fc in Eq. (fTSi) . 
where v denotes the density of states (assumed constant), we obtain the well-known 
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Meir-Wingreen expression for the current [I56I57II 

TiT. 



Here, 



I =2e ^ / dekAd{ek] 
i 1 + i _i J 



(23) 



A,{e,) = —J2ln^[GT,A^k) 



(24) 



is the nonequilibrium spectral function of the dot and is in general a function of 
the bias voltage $. However, in our notation we shall suppress this explicit bias 
dependence of the spectral function for convenience. In deriving Eq. (124)) we have 
adopted the standard procedure and linearized the spectrum about the Fermi sur- 
face. Furthermore, Fq, = Trt^z/ denotes the partial broadening of the level due to the 
coupling to the lead a. 

Next, we derive a more general form of the Meir Wingreen formula. This is accom- 
plished by the following sequence of steps. 



Proceeding as in Eq. (f2T|) one observes 



t [iU^n 



t-^ Ulr, 



OnO+ 



iUn — C + Cy 

-d. 



(25) 



Using a straightforward manipulation on the lines of Eq. (1221) . the first part of the 
expression in Eq. (l25l) can be simplified to give 



.0+ 



C + Cy 



(26) 



Thus, using Eq. (1181) together with Eqs. (1221) and (1261) we get 



I = ^ya^(g. t (T = 0)-g .t(r = 0)) 

/5 \ Ufc + ^ - ZWn - £ + £y J / J 
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Using an identity from Appendix IC.ll and integrating over in the last 2 terms, 
we find 



oo / <T) - 



(28) 



Now, 2z7r^,„,t(r = 0) = 2m (4d.) = 1'kG< ^,{t = 0) = rfe,G< (e,). 



Thus, we finally recover the more general form of the Meir-Wingreen formula: 



/ = y] otVa_ f deklf (ek - a— 
27r J-Qo [ V 2 



(29) 



We proceed in a similar fashion to derive an expression for the charge occupation 
on the dot (See Appendix C for details) 



where 



1 



f ek 



^ J dekAd{ek) 
J defcA,(efc)r«'(efc,$) 



f(ek + 



(30) 



(31) 



Note that on taking the limit $ — )> 0, Eq. (l30l) reduces to its well-known form 
from equilibrium many-body theory. In Appendix IC.ll we provide a straightforward 
derivation of this result using the spectral representation within the effective equi- 
librium formulation. This result can also be obtained within the Keldysh framework 
as we have shown in Appendix IC.2[ 

We thus identify the non-equilibrium spectral function Ad{uj) as the central quan- 
tity of interest, which one can use to calculate transport observables. In the effective 
equilibrium scheme this can be done directly, without resorting to the Keldysh con- 
tour which leads to Dyson equations which couple the different Green's functions. 
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3. 3 Simple Application: Non-Interacting Resonant Level 



The situation in which interactions on the dot are absent {i.e., Hmt = 0) corresponds 
to the double-barrier problem in quantum mechanics and presents the simplest sce- 
nario, and has been solved exactly via numerous approaches I162I63L We shall re- 
call the key results for this model outlining the methodology we adopt for obtaining 
them within the effective equilibrium framework. In Section 4 we shall use these 
results as the starting point for perturbative calculations for interacting theories. 

The formal expression for the Lippmann-Schwinger [I59I60II operators is given by 



1 



(32) 

,(0) 



Cafc - + 17] 

For the case of the non-interacting resonant level we shall use the notation ip^ka 
to denote the Lippmann-Schwinger operators in order to distinguish them from the 
Lippmann-Schwinger operators of interacting models. It is straightforward to show 
that 



,/,(o)t 



9d{ek)dl + -^2^ 



follows immediately from Eq. (fTTI) . Here, we have defined 

gdi^k) 



^a'k'a 



(33) 



(34) 



efc -ed + iV 

and r = Fi + F_i denotes the electronic level broadening of the dot. 

Eq. (|33l) can be readily inverted to express {c} and {d} in terms of {ij)^^^} to give 



t 



Y.9:iek)^^l- 



ak 



(35) 



cxkcT raka q / ^ 



(0)t 



9d{^k' 



,/,(o)t 

ra'k'c 



In the imaginary-time formalism we use Eq. (|23l) to compute the current 

et 1 



Im 



akuLjJri 



Using Eqs. (|35l) and (|36l) . the current can be finally written as 



et^ 1 
/ = Im 

n (3 



J2 J2^9d{ek')GAo) Ao)^{iuJn) 



akauJn Oi'k' 



(36) 



(37) 



(38) 
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Note that the expression for the current, in terms of the imaginary-time Green's 
function Q (o) .m {iun) is exact even when including interactions on the dot. In 

^ a' k' a ^otka 

this context however, the density matrix is no longer diagonal in terms of {^/'*^°^}, 
the Lippmann-Schwinger states of the non-interacting resonant level. Instead, they 
are diagonalized by the Lippmann-Schwinger operators of the interacting model 

The charge occupation on the dot is given by 



For the non-interacting resonant level the Green's functions in terms of Lippmann- 
Schwinger operators are trivial to evaluate, since the action in the generating func- 
tional is quadratic and diagonal in terms of these states and hence 



.0+ 



^ a'k' <t' ^aktr 



«2 



(40) 



Using this expression in Eq. (1381) gives 

foo 1 



TT J-oo (e - ed) 



Y2 



de. 



(41) 



Comparison of Eq. (|4T]) with Eq. (l23l) allows one to extract the nonequilibrium 
spectral function 



AAe) 



2T 



(42) 



It is instructive to note that for the non-interacting case the spectral function is inde- 
pendent of the bias. The introduction of interactions generally imparts to the spec- 
tral function a nontrivial bias dependence, as exemplified by the Anderson model 
which we shall discuss in Section 5. 



The electron occupation on the dot is evaluated in a similar fashion, namely 

X /"OO 



nd 



2 J — oo 



AAe] 



f{e 
1 



2 



TT J-oo (e - EdY + 



2" 



de 



f e 



de. 



(43) 



In the noninteracting case it is possible to express the bias operator Y explicitly in 
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terms of the electron operators in the leads and dot, 



akcr 




( 



t 




{gd{ek)dicaka + h.c.) 



Caka + h.C. 



) 



(44) 



which is far more intricate than the chemical potential term in the initial density 
matrix Yq = J2aka ^c^aka'^aka, sincc it encodes the entire nonequilibrium bound- 
ary condition of the system, and the bias voltage manifests itself as off-diagonal 
contributions for the electrons in the leads and the dot. This can be attributed to the 
highly non-local nature of the scattering states. 



4 Perturbation Theory for Interacting Models 



In this Section we develop a novel perturbative framework for interacting theories 
which facilitates the use of well-established many body techniques. 

The Section is structured as follows. In Subsection 4.1 we obtain a formal expan- 
sion of the Lippmann-Schwinger operators of the interacting system in terms 
of the Lippmann-Schwinger operators of the non-interacting resonant level {'0^°^}, 
in powers of the interaction strength. This expansion is used in Subsection 4.2 to ex- 
plicitly construct the ejfective Hamihonian, and hence the density operator. In Sub- 
section 4.3 we lay the framework for computing real and imaginary-time Green's 
functions, perturbatively in the strength of the interaction, using the effective equi- 
librium density operator of Subsection 4.2. 

A key feature of this formalism is that it is an effective equilibrium theory. It by- 
passes the need to define Green's functions on the Keldysh contour and eliminates 
the necessity to invoke coupled Dyson's equations. However, in contrast to the 
Keldysh formalism, the effective Hamihonian is not available a priori and needs 
to be explicitly constructed to the desired order. 

It is important to note that the scheme presented in this Section is a model inde- 
pendent generic procedure, barring the fact that the interactions are assumed to 
be confined to the vicinity of dot. In Section 5, we use the Anderson model as a 
prototype to illustrate the details of this scheme. 
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4. 1 Construction of Lippmann-Schwinger operators 



In this Subsection we sketch the formal expansion of the Lippmann- Schwinger 
operators, in powers of the interaction strength, using a convenient starting point. 
Typically the point of reference is an exactly solvable model, for example the non- 
interacting level, the Toulouse point of the Kondo model or the infinite-U Anderson 
model in the iV — oo limit {N being the number of "flavors" of the electrons). For 
the purpose of the ensuing discussion we use the non-interacting resonant level as 
our point of reference and construct the Lippmann- Schwinger operators in terms 
of {ip^'^^} in powers of the interaction strength. 

We introduce the additional superoperators C and Cx such that their action on any 
operator O is given by CO = [H — Hiau O] and CjO = [Hi^u O] respectively. 
Using Eq. (fTTI) we obtain 



-'akc 



+ 



1 



ek- C + irj 



-"aku 



+ 



1 - 



1 



Cfc - C + 17] 



/=0 L 



ek - C + 17] 



ek - C + 17] 
1 



ek - C + 17] 



^TCaka- 



(45) 



Since the Lippmann-Schwinger operators for the non-interacting resonant level sat- 
isfy 



aka ah 

we can rewrite Eq. (|45l) as 



(o)t _ J 1 



ek - C + 17] 



aka ' 



(46) 



i'^aka — "^akt + X! 



1=1 I 



ek - C + i7] 



Waka ^aka)- 



(47) 



This expression can be simplified to give 



i'aka — i'akl + X] 



2 oo 



1=1 



ek - C + i7] 



E (48) 



J a'k 



ek - ek' + ir] 



This can be expressed symbolically as 



i'aka — E '^aka 



(49) 



where the superscript / denotes the contribution proportional to the l-th power of 
Hint- 



is 



The formal expansion in Eq. (|49] ) is used in the following subsection to construct 
the density matrix order by order in the interaction strength. 



4. 2 Construction of the Density Matrix 



The steady state density matrix given by Eq. ([8]) is of the form 



P 



-f}{H-Y) 



(50) 



where 1-L = H — Y denotes the effective Hamiltonian for the system. The effective 
Hamiltonian is diagonal in terms of the Lippmann-Schwinger operators and can be 
written as 



(51) 



Our objective is to construct the effective Hamiltonian, and thereby the density ma- 
trix, in powers of the interaction strength. This will be instrumental for computing 
observables via nonequilibrium Green's functions, as we shall illustrate in Subsec- 
tion 4.3. Inserting Eq. ( |49l ) in Eq. (|5T1) above, we get 



oka ni,ri2=0 



E 

m=0 



Eh-«^)E^&ir^^ 

aka ^ ^ ^ p=0 



m=0 



(52) 



Here, 



akcj ^ ^ ^ p=0 



(53) 



denotes the contribution to the effective Hamiltonian H proportional to (iJir 
Note that the m = term. 



n^'' = E ( 

aka 



^ 2 J ^Cffccr^afca 



(54) 



is precisely the effective Hamiltonian for the non-interacting resonant level. 



We emphasize the fact that the perturbative scheme developed here is very distinct 
from the usual perturbative expansion familiar from equilibrium quantum field the- 
ory. In the present approach, the "perturbation" is not available a priori, rather 
the effective Hamiltonian is to be computed to the desired order in the interaction 
strength before resorting to diagrammatic methods. 
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4. 3 Perturbative computation of Green 's functions 



Next, we use the effective Hamiltonian from Eq. (|52l) to compute Green's functions 
to any desired order in the interaction strength. We present two distinct frameworks 
for computations, in real and imaginary-time respectively. 

The imaginary-time formulation involves a straightforward implementation of the 
effective Hamiltonian via a functional integral representation. This is achieved 
by a coherent state representation (over Grassman variables) of the Lippmann- 
Schwinger operators of the non-interacting resonant level. The motivation behind 
the use of this representation is the fact that these operators are simply related to 
the {c} and {d} operators, as given by Eqs. (l35l) and (|36l) . All transport quantities 
are formulated in terms of the {c} and {d} operators which represent the original 
degrees of freedom of the system. 

For the {?/;*^°^} operators, the generic single-particle Green's function is given by 

^^(0) ^(o)t (^) = -{t [<k., (r) V^SL, (0)] ) . 

Note, the density matrix used in computing the above expectation value is diagonal 
in terms of the Lippmann-Schwinger states of the interacting model. We apply 
the formal expansion for H. from Eq. (l53l) in the above expression, in order to 
rewrite the effective Hamiltonian in terms of the Lippmann-Schwinger states of 
the non-interacting resonant level retaining terms to the desired order in Hi^i- This 
is implemented by expressing the partition function as a coherent state functional 
integral 



Q^w ^(o)t (r) 



1 



IX 



exp 



a.ka 

- f dr' E + -H (^(°)*, 



ako 



(55) 



Here, the partition function is given by 



Z =tr 

X exp 



ako 



Oikcr 



(56) 



and we have suppressed explicit time labels for the Grassmann variables featuring 
in the effective Hamiltonian. Using Eq. (l53l) in the exponent of the above equation 
we can formally write 
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,,(0)t {t ) 



1 



X exp 



^ aka I m=0 



(57) 



Introducing the free action S' = /^f dr' Eaka {^it (^')^r'^il(r') + ^^^^ (v^®*, ^(°)) } 



and using the abbreviated notation T-L^"^^ 
the exponent to obtain 



(0)* 



'H^'^)[t'), we can expand 



1 



i/=0 i=l''^ \m=l 



fji:^ H (/^r,««™.'(r,))...(/''dr,«'"-.)M 

n=0 I «=1 mi+...+m,=n / V"'^ / 



(58) 



In the last step, we have reorganized powers of i^int by introducing the indices 
mi, 7712, ■ ■ ■ ,mi G {1,2,...} and used the convention that the term in curly brackets 
is unity for n = 0. This allows us to isolate the order n contribution. Thus we can 
write 

oo 

V\ ^(o)t W = E^i(o) ^(o)t (r), (59) 

f-Qifeio-i f-o:2fc2<^2 ^aifcio"! V'a2*2°"2 

where the n-th order contribution is explicitly given by 

,(o)t (r)= ^/n(^^€i:^^€2j€2L2(o)<Uwe-'^ 



1=1 ( -ini+...+mi=n 



The formal expression in Eq. (l60l) above, captures the contribution proportional 
to (i/int)" to the single-particle Green's function. This procedure can be used to 
compute Green's functions, and hence observables, perturbatively. The imaginary- 
time formulation enables one to use the gamut of functional methods and may also 
be used for nonperturbative analyses, e.g. the infinite-t/ Anderson model in the 
limit of — !■ oo. 
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We now illustrate the procedure for computing Green's functions in real-time. In 
this case the expansion is more intricate, since it involves a double expansion in 
powers of Hint- In Section 3.2 we established the retarded electron Green's function 
of the dot as the quantity of central importance in the real-time description. This 
can be expressed in the frequency domain as 



1 



uj + C + ir] 



tr 



tr [e-/3«] 



(61) 



Note that the bias operator appears in both the exponent and the Liouvillian super- 
operator. To compute the above expression perturbatively, we systematically collect 
powers of U arising from the expansion of both the exponent and also the Liouvil- 
lian. First, we expand ^^l^j^^^ d^ in powers of U 



DJuj) 



1 



u + C + ir] 



ni=0 



ni=0 

which implies the recursion relation 



ni 



UJ + C + ir] 



oj + C + ir] 



u + C + ir] 



dn 



(62) 



(63) 



For any observable O the expectation value can be expanded in powers of the in- 
teraction as 



tr 



-mo 



OO 

n2=0 



i=l 



, m=l 



"2 



Er 



ir n/„ *. 



X Yl HI"">(ti)...H<"'">(t,) 

\ mi+...+mi=n2 



E(^) 

n2=0 



("2) 



(64) 



where mi, m2, . . . ,mi E {1,2,...}. Here, (C)"^ represents the n2-th order con- 
tribution to the expectation value of the observable. We use the convention given 
below Eq. (l58l) for the sums. The subscript "conn" implies we discard all "discon- 
nected" diagrams. Once the perturbation expansion is systematically implemented 
we use the standard diagrammatic rules for computing the contributions of the in- 
dividual terms of the expansion. Finally, we systematically collect powers of Hint 
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arising from the two distinct expansions 



ni,n2=0 



We rearrange the double sum and group together terms containing identical powers 
of Hint to obtain 



7i 



n=0 lk=0 

oo 

n=0 



(n-fe) 



(66) 



■yret(n) 

dad^a 

function of the dot. 



Here, G^'^ } denotes the n-th order contribution to the retarded electron Green's 

dad' 



We have thus constructed an explicit perturbative expansion for imaginary-time 
and real-time Green's functions applicable to generic interacting quantum impurity 
models. When compared to the Sch winger- Keldysh approach, this scheme allows 
a direct computation of the Green's function and bypasses the need for the Keldysh 
contour. 



5 Anderson Impurity Model 



In this Section, we use the Anderson impurity model as a prototype to implement 
the perturbative framework presented in Section 4. The Anderson impurity model 
||54| describes a commonly encountered experimental scenario and serves as the 
microscopic model for a wide variety of fundamental physical phenomena. The 
interaction is given by 

^^int = (firf - 1) = ^ ^ d\d}_„d-„d„, (67) 

and denotes the Coulomb interaction between electrons of opposite spin projections 
on the quantum dot. Here, rifi = Y^a dl.d^ denotes the electron number operator of 
the dot. Thus, the interaction can be rewritten as 

^ V^V 1,2,3,4,(7 

where we have introduced the abbreviated notation / = [oLiki] and gi = gd{eki), 
implying a sum over the lead index a and the quantum number k. Furthermore, we 
define Y.i = Y^aiki for the subsequent discussion. 
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Our analysis of the Anderson impurity model is organized as follows. We begin 
with a brief discussion of the mean field results of the Anderson model which fol- 
lows immediately from the results obtained for the noninteracting resonant level in 
Section 3.3. In Subsection 5.2 we calculate the self energy to first order in the inter- 
action strength within the imaginary-time framework and compare it with the mean 
field results. In Subsection 5.3 we study the nonequilibrium dynamics of the An- 
derson impurity at the particle-hole symmetric point within the second-order Bom 
approximation. We also consider deviations from the particle-hole symmetric point 
by introducing a local magnetic field on the dot. 



5. 1 Mean Field Theory 



As a preliminary analysis we treat the Anderson model out of equilibrium within 
the mean-field approximation. This amounts to replacing the interaction Hamilto- 
nian above with 

H^^ = UY.{n^Md.. (69) 

(T 

Here, n^j = d)^d^ gives the electron number operator of the dot for spin projec- 
tion a. At the mean-field level, the interaction renormalizes the dot energy — )■ 
[ed + U ifi^a)), such that 

^r = E(e^ + f^(^-))4rf- (70) 



Using Eq. (|4T1) . the expression for the current reads 



27r v^-oo {e-td-U{h^,)y + r^ 



rfe, (71) 



where the occupation number of the dot is determined self-consistently II641I via 



27r 



1 



(e-e,-?7(n_))2 + P 



2" 



de. (72) 



One observes that the condition for particle-hole symmetry, i.e., {n„) = 1/2 is 
realized for = — 1//2, as confirmed by first order in perturbation theory. 



21 



5.2 First-order analysis 



The first order contribution to the Lippmann-Schwinger operator, cf. Eq. (|49l) . is 
given by 



^S] = Ug,{e,) ( - ) E 



919293 



n 



123 



efc - €1-62 + 63 + if] 



^'S'^flUfl. (73) 



Using this expression in conjunction with Eg. (1531) we construct the effective Hamil- 
tonian to first-order in the interaction 



aka 



(o)t„/,{i) 



aka 



(74) 



121'2'cr 



where we have defined the matrix element 



;i2|V| 1'2') = u 



ev - ei - {ay - ai) 



9i929v92' 



ey + £2' - Ci - 62 + ZT] 



. (75) 



Note that (12 |V| 1'2') is not symmetric under (1 ^ 2, T ^ 2'). In diagrammatic 
language, one needs to distinguish opposite ends of the (two-body) interaction line, 
as we have illustrated in Fig. [21 



Assuming identical leads, we use Eq. (1381) to compute the first order correction to 
the current (using standard diagrammatic techniques). 



/(i)=_^llm 



Vt (3 



Un 



(0) 



TT 



dei- 



ei-e,) /(ei-f -/(ei + f 



[(ei-e,)2 + P] [(ei-e,)2 + p 



(76) 



where 

denotes the charge occupation of the dot for the noninteracting resonant level as 
given by Eq. (l43l) . 

Note, for the first order computation only the Hartree diagrams contribute whereas 
the Fock diagrams are irrelevant. Comparing the expression for the current in Eq. 
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(|76|) with Eq. (1231) we extract the spectral function to first order in the interaction 
2T 

A,{e) = - 

TX 

One observes that this result agrees precisely with the mean field solution discussed 
above, to first-order in U, in accordance with Hartree-Fock theory. 



5.3 Dynamics about particle-hole symmetric point 

Next, we study quantum transport through an Anderson impurity which is at the 
particle-hole symmetric point, given by = — ^ . The special feature of the particle- 
hole symmetric point is that by tuning the interaction strength U, one can pass 
smoothly from the weak to the strong coupling regime such that the Fermi liq- 
uid fixed point of the Hamiltonian remains invariant. Also, the Anderson impurity 
model maps exactly onto the Kondo model in the limit U — )■ oo at this point in 
parameter space. 

In the following discussion it is convenient to absorb the bare energy of the dot 
in the interaction term, i.e., we let — )■ q = + ^, thereby generating an ad- 
ditional term — ^ d^da- The particle-hole symmetric point is thus specified by 
€d = 0. This redefinition assures that the pole of the bare propagator is correctly 
positioned, and becomes particularly important when we extend our results to the 
strong coupling regime. The interaction term is consequently redefined as 

Hint = ^i^d- if ■ (79) 



1 

[(e-e,)2-FP 



(78) 



Only even powers of H-^at feature in the perturbative expansion of the spectral func- 
tion at the particle-hole symmetric point, and the leading correction is of order 
O (f/^). In our discussion, we adopt the real-time approach to compute the spectral 
function Ad{uj). In this context the real-time approach is computationally simpler 
than the imaginary-time approach since the n-th order calculation in U requires us 
to compute the effective Hamiltonian only to order n — 1, at the cost of evaluating 
L'i")(tc;),cf.Eq. 

The non-interacting part of the Green's function is identified to be 

which immediately followed from Eq. (|66l) . Evaluating the second order contribu- 



23 



tion we obtain (see Appendix iDl for details) 
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w + ei - 62 - 63 + z?7 
where we have abbreviated = i / (^n — f ) + / (^n + f ) 



(81) 



The Dyson equation for the retarded electron Green's function of the dot is given 
by 



Tet(0)| 

dcrdt 



(82) 



Here, Y7f{uj) is the proper self-energy obtained exclusively from contributions of 
one-particle irreducible diagrams. The Dyson equation above can be rewritten in 
the alternative form 



ret(O) , 



(83) 



where T.^^\uj) represents the complete self-energy of the system, obtained by sum- 
ming the 1 -particle irreducible (IPI) diagrams. 

Expanding G^f ^ (w) and SJ^'(a;) in powers of iJint it is possible to relate the n-th 
order {n>l) contribution to the self-energy 



dad' 



dad' 



(84) 



Since the first order result is zero, we conclude that (2) = (2). The second- 
order contribution to the self-energy is directly read off from Eq. (ISTI) . and can be 
simplified to give 



^ret(2) I 
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u + 3iT 
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(85) 



The low-bias behavior for the entire range of the parameter U is captured in the 
Fermi liquid expansion ll55l for the non-equilibrium spectral function and the dif- 
ferential conductance. This is typically expressed in terms of the susceptibilities 
Xn and Xtt as [1551 



Ad{uj) = — 
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(86) 



24 



and 



G — G(\ 



V 



V 



(87) 



Here, G = ej^ denotes the differential conductance, and Gq = 2e'^/h is the con- 
ductance quantum. 



The susceptibilities xu Xtt analytic functions of U/ (vrF) and their analytic 
properties have been extensively studied. It is well-known that their Taylor expan- 
sions display a rapid convergence |[65l . Extracting the spectral function from the 
expression for the Green's function to second-order in Eq. (|85l) . one obtains 
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7rf 



and similarly 
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l/Uy/irTy 3 
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(88) 



= Go 




(89) 



This agrees with Eqs. (l86l) and (l87l) where we use the well-known perturbative ex- 
pansions for the susceptibilities = (^) andxtt = 1+i (6 — ^) 




Strictly speaking, the expression for the self-energy to order f/^ is a well-controlled 
approximation to the actual self-energy only when U <^ T. On the other hand 
Kondo physics becomes relevant in the regime U ^ T. For a qualitative analysis 
of this regime we apply the Bom approximation using the second-order self energy 
obtained above ll66l . 

It is in general difficult to rigorously justify the exact bounds of validity of the Bom 
approximation without identifying an explicit control parameter, but this scheme is 
seen to give physical results for fairly large values of U (~ SF) at the particle-hole 
symmetric point [|66ll . 

When U ^ T and $ ^ F, we notice the appearance of the Abrikosov-Suhl reso- 
nance centered at w = reminiscent of the Kondo resonance, which corresponds 
to the case U — oo. This many-body resonance needs to be distinguished from 
the sidebands corresponding to the bare atomic levels which continue to persist. 
The resonance at the Fermi level diminishes with an increase in the bias voltage 
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as shown in Fig. 3. One may conjecture that the effect of finite voltage might be 
captured by an effective temperature; this is in some sense related to a dephasing 
rate induced by a finite current [|68l . Our results corroborate the absence of a two- 
channel Kondo effect that would be associated with the coupling to two fermionic 
continua at different chemical potentials II69I70II . At large voltages the resonance 
is completely destroyed and one obtains the residual atomic sidebands at ± 
Comparing Fig. 3 and Fig. 4 we also confirm that the central peak has a diminishing 
width when increasing U . 

This second-order result, which is exact in the bias voltage, reduces in the low-bias 
and high-bias limits to those in similar discussion using the Schwinger-Keldysh 
scheme II66I67I71II . Qualitatively they agree with the NRG results and are con- 
firmed by recent experiments [I22I49II . The splitting of the Kondo resonance at in- 
termediate biases is still an open issue since other groups have observed a split in 
the Abrikosov-Suhl resonance with increasing bias using a fourth-order calculation 
||72l . However, recent numerical calculations in the intermediate coupling regime 
[I22I37I39I42I58I731 indicate that the peak splitting from the fourth-order perturba- 
tion theory has been overestimated. 

Finally, we study deviations of the system from the particle-hole symmetric point. 
For this purpose we introduce a local magnetic field {H) on the dot which implies 
— ^ + thereby breaking particle-hole symmetry at q = 0. Here, the devi- 
ation 5a = — (Ty is given in the units defined hy g fib = 1, where g is the electron 
g-factor and /i^ denotes the Bohr magneton. The retarded electron Green's function 
of the dot for the non-interacting level now differentiates between the different spin 
projections 

'^'^"■'^ U — Oct + tl. 



Away from the particle-hole symmetric point certain diagrammatic contributions to 
the self-energy which earlier evaluated to zero, cease to do so. The leading order 
contribution to the self-energy is at first-order. Using Eq. (|66l ) we have computed 
the self-energy to second-order in the interaction, with details summarized in Ap- 
pendix |El The expression for the self-energy correctly captures the entire bias and 
field dependence to a given order in U. Using this result to approximate the retarded 
electron Green's function of the dot we observe a split in the Abrikosov-Suhl res- 
onance as we increase the field strength H as shown in Fig. 5. This split becomes 
discernible when i/ ~ F. Furthermore, the sidebands shift further apart as ex- 
pected. For the Kondo model, which corresponds to the extreme limit U ^ oo, the 
effect of a magnetic field has been rigorously studied II74I75I76II . 

The current in Fig. 6 is obtained via Eq. (l23l) . which at T = corresponds to 
integrating the spectral function over a window of length $ centered on the dot. 
Note that the spectral function is a function of the bias and we observe suppression 
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of the Abriksov-Suhl resonance with increasing bias accompanied by the atomic 
sidebands shifting further apart. In the $ — )• oo limit the spectral function is in- 
tegrated over the entire range of energies, which by the sum rule / duiAd{ui) = 1 
assures that the current saturates to the same asymptotic value. We see suppres- 
sion of the current with increasing U for low biases corresponding to a sharpening 
of the Abriksov-Suhl resonance. The sudden increase of the current for higher $ 
corresponds to the contribution of the sidebands. 



6 Conclusion 

In this paper, first we have discussed the theoretical foundations of the effective 
equilibrium theory initiated by Hershfield ll50l . and demonstrated a correspondence 
with the observables computed using the Schwinger-Keldysh formalism. We have 
pedagogically re-derived the effective equilibrium density matrix approach pro- 
posed by Hershfield using the concept of an "open system", which does not require 
inclusion of additional relaxation mechanisms. Furthermore, we have illustrated 
the methodology of this formulation for a quantum dot, in which the reservoir 
leads play the role of (infinite) thermal baths, such that a unique steady state ex- 
ists. The latter can be cast into the form of an effective equilibrium density matrix, 
where the associated modified Hamiltonian can be explicitly written in terms of the 
Lippmann-Schwinger scattering state operators. 

We also introduced a systematic perturbative framework for interacting systems 
out of equilibrium within the effective equilibrium theory. This procedure involves 
a novel expansion in terms of the Lippmann-Schwinger states of the system, in a 
manner very distinct from the Schwinger-Keldysh formalism. Using this method 
the need for the Keldysh contour is eliminated and the Green's function can be 
computed directly. This allows use of the usual Feynman diagrams albeit with a 
modified Hamiltonian. Furthermore, the perturbation theory is free of unphysical 
infrared divergences. However, the structure of the effective interactions becomes 
more complex with increasing order in perturbation theory. 

We have illustrated our scheme with an application to the Anderson impurity sub- 
ject to a local magnetic field and/or a finite bias. The non-equilibrium electron spec- 
tral function is evaluated analytically to second-order in the interaction strength, 
while capturing the entire bias and magnetic field dependence. We use the sec- 
ond order Born approximation for the self-energy to simulate the strong-coupling 
(Kondo) regime. The fate of the Abrikosov-Suhl resonance as a function of the 
bias and magnetic field was inferred from this result. We observed that increasing 
the bias causes a gradual decrease in the resonance until ultimately only the two 
atomic sidebands remain, indicating that the Kondo peak is destroyed with bias 
voltage without a perceptible splitting of the resonance peak. This result reduces in 
the low-bias and high-bias limits to that in similar discussion using the Schwinger- 
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Keldysh scheme II66I67I71I I. Applying a magnetic field however causes the peak 
of the resonance to split, while the sidebands get farther apart due to the fact the 
singly occupied level is lowered in energy by —H/2 while the energy of the dou- 
bly occupied level increases by H/2. This agrees qualitatively with numerical and 
experimental results [122137 139I42I49II . As a final check of this method, our results 
reduce to the well-known Fermi liquid expansion at low bias, energies, temperature 
and magnetic field ll55l . 

The scheme has a potential for being generalized to the non-perturbative regime us- 
ing well-controlled techniques such as the large-N expansion [I34I77L and such ex- 
tensions are currently under investigation. In this case the Lippmann-Schwinger op- 
erators we use as our point of reference follow from an effectively non-interacting 
Hamiltonian, which arises from the large- saddle point solution of the functional 
integral obtained using the slave-boson formalism. There is scope to treat more so- 
phisticated systems such as Luttinger liquids within this formalism by combining 
bosonization and scattering techniques. Also, a possible extension of the formalism 
to include thermoelectric transport in addition to bias-induced transport, is being 
explored. We emphasize that the effective equilibrium approach is that it is con- 
ducive to the implementation of numerical methods, such as numerical RG I122I23II 
and Quantum Monte Carlo [|42]| . 

This method is independent of the exact nature of the interaction, provided it is 
localized in the vicinity of the quantum impurity, and can be used to treat other 
models such as the interacting resonant level. Furthermore, the extension to multi- 
ple leads and other more complicated geometries can be straightforwardly accom- 
plished. 

Acknowledgments: We acknowledge N. Andrei, H. U. Baranger and T. L. Schmidt 
for stimulating discussions. This work is supported by the Department of Energy 
under grant DE-FG02-08ER46541 (RD. and K.L.H.), by the Yale Center for Quan- 
tum Information Physics (NSF DMR-0653377, J.K. and K.L.H.) and by the NSF 
under grant DMR-0907150 (J.E.H.). 



A The Y operator 

In this Appendix we present an alternative derivation of the steady state density 
matrix. To build the F-operator, the starting point is the expansion of the steady- 
state density matrix into a power series in the tunneling Hamiltonian Ht- This can 
be accomplished by employing the time-dependent framework of the open-system 
limit where the steady-state density matrix is obtained by adiabatically switching 
on the coupling in the far past. One obtains analytical expressions for the power 
series in Ht- Order for order comparison of these expansions then leads to a system 
of nested differential equations for the bias operators. Finally, the solution to this 
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system of differential equations is then shown to be identical with the representation 
of Y in terms of Lippmann-Schwinger operators, see Eq. (fT2l) . 

We should emphasize the fact that the expansion in powers of Ht is a purely formal 
procedure, which we adopt for the sake of a systematic comparison, and the proof 
below is non-perturbative in Ht- 



A.l Expansion of p using the effective equilibrium representation 



In the derivation of the explicit form of the Y operator, it is convenient to collect 
terms in orders of the tunneling Hamiltonian Ht- We thus start by expanding Y = 
Y.'^=o^n into a series in powers of the tunneling, such that F„ oc (Ht)"'. In the 
following, the index n will always be used for power counting of the tunneling 
Hamiltonian Ht- From Eq. ([8]) we thus obtain 



p = exp 
= exp 



-/3 |(i7o - Yo) + (Ht - - ^ F„ 

oo 



n=2 



n=0 



(A.1) 



Here, we have regrouped the Hamiltonian with the Y operator order by order in the 
auxiliary operator Xn ~ {Ht)"', which is hence defined as 



ifo - >o, n = 
X„ = Iht-Yi, n = l 
-Yn, n > 2 . 



(A.2) 



One can expand the exponential operator to collect terms in powers of the tunneling 
Ht such that 



p = E 



I oo 



1=0 n,...,i;=0 

oo oo /' Q\l 



EE 

n=0 1=0 



E -^ii • • - ^i, = E 



(A.3) 



iiH \~ii=n 



n=0 



where, following our general notation, p„ denotes the order {HtY contribution to 
the steady-state density matrix. 
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A.2 Expansion of p using the open-system approach 



Let us now employ the open-system approach (as outlined in Section 2). In this 
case, the steady-state density matrix is obtained by switching on the tunneling in 
the early past to < 0. Up to this time t^, the leads are decoupled from the dot. 
The adiabatic switch-on of tunneling is facilitated by using Hxe^'^O^t — to) for the 
tunneling Hamiltonian (note that in this scenario, the total Hamiltonian is therefore 
time dependent). At time t = 0, transients have decayed and time evolution has 
turned the original density matrix po = p(t = to) into the steady-state density 
matrix p = p(t = 0). We note that the condition of adiabaticity is strictly true only 
in the limit rr-r ^ ?7, where Itol — ^ oo, which is assumed in the open-system limit 
11181 . We emphasize that the actual evaluation of this limit is deferred until the very 
end of all calculations. Keeping l/|to| and t] small but nonzero in the interim is 
crucial for mathematical clarity and for avoiding ill-defined expressions. 

This time, it will be convenient to work in the interaction picture (with respect to 
the tunneling), where in general 

Oi{t) = e^^o^^-^f^Oe"'^"^*"*"). (A.4) 

denotes the interaction picture of the operator O. Note that the time to (not t = 0) 
has been chosen as the reference time where Schrodinger and interaction pictures 
agree, O = C/(to). [In this we differ from the conventions adopted by Hershfield 
[|50l who chose different reference times for the Heisenberg and interaction repre- 
sentation.] 

In the interaction picture, the density matrix satisfies the evolution equation 

tj^pi{t) = [HTAt),Pi{t)]. (A.5) 

This is formally solved by 



n=o k i=l 

oo 



rt 

dti 

'to 



[HTAtl),[HTAt2),[. . .[HTAtn),Po]] 



J2PnAt), (A.6) 



n=0 



where the = term simply fixes the boundary condition p(t = to) = po and T 
denotes time-ordering of operators. 

It should be noted that Eq. (IA.6I) already has the form of a power series in the 
tunneling and thus has been used to define the n-th order contribution p„,/(t) to 
the interaction-picture density matrix. For later purposes it is useful to note that the 
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contributions can alternatively be obtained from 

i.e., a system of nested differential equations associated with the boundary condi- 
tions po,/(^o) = Po and Pn,i{to) = for n > 1. 



A. 3 Derivation of differential equations for Yn 

For the results to be consistent, we require (assuming steady state is reached Il53l ) 
that the steady-state density matrix p be identical to the steady-state density matrix 
obtained in the open-system limit. This allows one to determine the correct form of 
the Y operator, order for order in the tunneling. 

We thus impose the identity of the steady-state density matrices 

where we have transformed p„ into the interaction picture. Eq. (IA.8I) is expected 
to hold for all t > 0, since at that point the interaction has been fully switched on, 
and the time-independent Hamiltonian (used in Hershfield's effective equilibrium 
approach) and the time-dependent Hamiltonian (from the open-system approach) 
are identical. 

We now differentiate the left-hand side of Eq. (IA.8I) with respect to time and use 
Eq. (IA.3D to obtain 

4„.w = f;^f^ E (A.9) 

1=1 '■• ii+-+i;=n k=l 

where Xi^^ j = ^^^oit-to) ^,^g-iHo[t-to) (jg^otes the interaction picture representa- 
tion of Xi^. Similarly, we may differentiate the right-hand side of Eq. (IA.8I) . The 
resulting commutator is given in Eq. (IA.7I ) and contains Pn-i,/(i)lZl for which we 
substitute the corresponding expression from Eq. (IA.3I ). This way, we obtain 



^ Note that due to the identity (IA.8I) we may, from here on, drop all bars on p. 
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J^Pnj{t) = i[pn-l,l{t),HTjit)] 



CO (_fj)l I 

^ E ^ E E x^At) ■ ■ ■ [x^,At)^ HrAt)] . . . x^t) 

1=1 ' 11+12+. ..+i;=n—l k=l 
oo (_g)l I 

^ E Sr^ E E ^nAt) ■ ■ ■ [x^^^iAt), HrAt)] ■ ■ ■ x^t), 



n 

=1 ■ 11+12 + . ..+j;=n k=l 



(A. 10) 



where for n < we define X„ = 0. In the last step of Eq. (lA.lOl) . the summation 
constraint is shifted from — 1 to to facilitate the comparison with Eq. (IA.9|) . 
This comparison yields the relation 



jXnj = i[X,,_^j{t),HT,i{t)]. (A.11) 

Finally, utilizing the relation (IA.2I) between X and Y operators, one finds that the 
Yn operators also satisfy 



jYr^j = 2[y„-i,/(t), HT,m- (A.12) 



To obtain the Y operator from these differential equations, it is crucial to specify 
the boundary conditions at the initial time t = to- Given the relation e'^''^*'^ = 
^iHo{t~to) ^~iHoit~to] ^ valid for any operator O, the boundary condition p(t = 
to) = Po implies [|53l 

limy,(t) = i^aiV„. (A.13) 

t\to Z „ 

It now remains to prove that the following interaction-picture expression of Y 



Yiit) =1 E »^ik.A^)^o.k.At)^ (A.14) 

akcr 



in terms of the Lippmann-Schwinger operators 'ipakaA'^) represents the solution to 
the above initial-value problem. Once we have recapitulated the crucial properties 
of these Lippmann-Schwinger operators in the following subsection, it will be sim- 
ple to confirm that this ansatz indeed solves the differential equation (IA.12I) subject 
to the boundary condition (IA.13I) . 
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A.4 Properties of the Lippmann-Schwinger operators 



It has been shown [[6TII for the generic model defined in Section 2 that the Lippmann- 
Schwinger operators are fermionic, 



and diagonalize the full Hamiltonian (including all interactions), 



(A.15) 



(A.16) 



ak 



One can expand ipak^. in powers of the tunneling Hamiltonian Ht, i.e., 

00 



(A. 17) 



n=0 



such that ip^akan ^ {HtT ■ In the interaction representation, the operators i^'^ako- 
satisfy 



d 
It 



subject to the boundary condition 



lim ^aka,l{t) = Caka- 
t\to 



(A. 18) 



(A.19) 



Eq. (|A.18I) can be further simplified, and collecting orders of {HtY', cast into the 
set of nested differential equations 

^^L.,n,/(^) = + Ai^ik.,n-iAt). HrAt)]- (A.20) 



As discussed in Ref. II6TI . the formal solution can be expressed compactly as 

1 



^1 



^ak J~' -r 01/ 



(A.21) 



From Eq. (fTTl) one can write the detailed form of the operator ipl^j^^ „ 

00 r / \ ™ "I 

V'L. = E 7^—-^^] ^E^kn- (A.22) 

n=0 IK'^ak - Lo + IT] J J „ 

Note that when the tunneling is set to zero, one finds indeed that ip'lf.^ — )■ cj^,^^. 
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A.5 Conclusion of the proof 



It is now simple to verify that the ansatz for Ynj{t) solves the set of differen- 
tial equations given by Eq. ( IA.14I) . and obeys the appropriate boundary condition 
(IA.13I) : We start by writing the (-f^r)" contribution to the Y operator as 



(A.23) 



Differentiating this with respect to time yields 



aka ^ p=0 L 



d 



-rAa,pAt) i^aka,n-pAt) 



dt 



d 



+ Vaka,pAt) \ ^i^aka,n-pAt) 



(A.24) 



so that in conjunction with Eq. (|A.18I) one obtains 



d $ 

aka p=0 

This concludes the proof. 

In summary, we find that the steady state dynamics of the system can be described 
by an effective equilibrium density matrix of the form 



where the operator 



p = e 



^ = Y E ai^ika'^ako 
aka 



(A.26) 



(A.27) 



encodes the entire nonequilibrium boundary condition of the system. 



A.6 Recursion relations for expectation values of observables 



In addition to the previous proof, we follow Hershfield [[501 and underpin the equiv- 
alence of the adiabatic approach and the effective equilibrium approach by showing 
that they lead to identical recursion relations for expectation values of observables. 
Again, the systematic use of the open-system limit makes the proof sound. 

Let O denote a generic observable. To derive the first recursion relation, we will 
decompose the steady-state expectation value (O) into a series counting the powers 
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of the tunneling Hamiltonian Ht- In the first step, we thus substitute the expansion 



(O) 



tl[pO] tljEZoPnO] 

tr[p] tr[E^=oP-]' 



(A.28) 



Pulling out a factor of 1/ tr[po]. expanding the denominator as a geometric series, 
and collecting terms order for order in Ht, we can rewrite this as 



oo n oo ; 

(0)=EE(-iy E' n 

n=OZ=0 ji,...,ji=l s=l 

oo 

= E(^)n, 

n=0 



tr[po 



tr[p, 



'"-El=i. 



O] 



(A.29) 



where Y.' denotes a restricted summation, subject to the condition I]i=i js < n. 
Eq. (IA.29I) allows one to prove the important recursion relation 



(A.30) 



which relates the n-th order term to the expectation value of O with respect to p„, 
and all lower-order terms {0)m (m = 0, . . . ,n — 1). 

Now, we turn to the time-dependent representation assuming an adiabatic switch- 
on of the tunneling. The steady-state expectation value of an observable may now 
be obtained via 

tr[po] 

where we have switched to the interaction picture. To arrive at the desired recursion 
relation, we again decompose this into a power series of the tunneling Hamiltonian, 



(A.32) 



where 



{0)r 



n! tr po 



tr rn 



i=l 



t=o 



dti 



to 



(A.33) 



X [HrAti), [HrAh), [• • • [i/r,/(tn), Po]] • • •]0/(O) • 



This expression for {0)n allows one to derive the second recursion relation, which 
reads 



{0)n 



tr[p„,j(O)0j 
tr[Po] 



(A.34) 



35 



The details of this derivation are given in Appendix|Bl The proof makes explicit use 
of the factorization of two-time correlation functions in the limit of large time sep- 
aration Il53l . The agreement between Eqs. (IA.30I) and ( IA.34I) underpins the equiva- 
lence between the time-dependent adiabatic approach and the effective equilibrium 
approach. 



B Derivation of the recursive form of O. 



In this Appendix we provide the detailed derivation of the recursion relation (IA.34I ) 
when working within the time-dependent approach, switching on the tunneling adi- 
abatically. Our starting point for the derivation is the expression for the order- (if^) 
contribution to the steady-state expectation value of some operator O, 

{0)n =^T^ tr { dtn dtn-l " " " / dt, 
tr[poJ [ J to J to J to 

[HrAtr,), [Ht,i{K-i), [. . . [ifT,/(ti),Po]] • • ]Oi{t = 0)|, (B.l) 

see Eq. (IA.33I) in the main text. Carrying out the ti integration and using Eq. (IA.5I) . 
one obtains 



tr[po] 



rtn rt3 

dtn / dtn-1 ■■■ dt2 (B.2) 

to J to •'to 



X [HtAQ, [. . . [i^T,/(t2), (piAh) - Pl,/(to))] . . -piio) 



We now retain the term involving pij{to), and proceed by carrying out the t2 in- 
tegration for the term containing pi,/(t2)- This integration yields 2 terms: one with 
P2,/(io) which we retain, and the other with p2,i{H) which we subject to further 
integration. Continuing in this fashion until all the variables ti, . . . , in the latter 
term have been integrated out, we arrive at the result 

[, ,(0)0,(0)] (3 3, 

tr[po] 

" i^-iY^P l'tn+i=0 rtn rtp+i 

— 2^ — r- / dtn, / dtn^i ■ ■ ■ dtp^i 

tr[poJ J to J to J to 

X tr |[i7T,/(tn), [. . . , [HTAtp+i),PpAto)] ■ ■ -PjiO)^. (B.4) 

Let us focus on the p-th term of the above sum. One observes that when p = n, 
none of the integrals are present and the argument of the trace is just pj, /(to)0/(O). 
Suppressing the integrals and other c-numbers we isolate the product of operators 
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within the trace 



ti ^[HtAQ, [..., [HTAtp+,),Ppj{to)]] . . .]O/(0)|. 

We can cycle Ppj{to) in the above expression out of the nested commutators, 
thereby inverting the nested commutator structure to re-express this as 



tr Ppj{to)[. . . [0/(0), HtAQ], ...HTj{tp+,)] . 



(B.5) 



This allows us to rewrite „ in the form 

tr [p„,/(0)O,( 



{Oh 



E 



tr[po] 



1 tr[p, 



'oj Jto 



dt 



in 

n I dtn—l 

to 



■p+2 



dt 



to 



p+1 



X tr 



To proceed further, let 



((••■))o 



tr [e-Pi.Ho-Yo) . . 
tr (e-/3(^^o-^o)) 



(B.6) 



(B.7) 



denote an expectation value with respect to the density matrix po- It has been shown 
by Doyon and Andrei in Ref. Il53l that, in the open system limit, long-time correla- 
tion functions factorize (this argument is rigorous at finite temperatures), i.e. 



mh)Oit,)))o ^ ((0(ti)))o((0(t2)))o 



(B.8) 



as \ti — ^2 1 — oo. We will now employ this factorization in Eq. (IB.6h . Renaming the 
integration variables t[ = tp+i, ^2 = tp+2, ■ ■ ■, t'n-p = and inserting poPo ^ = 1 
we get 



{0)n 

n 

-E 



tr [p„,,(O)0,(O)] 



(B.9) 



1 tr[pi 



'0 J to 



dt 



n—p 



dt' 



to 



n—p— I 



to 



c/t'itr<jpo Po Vp,/(^o) 

A(to) 



x[[...[Om,HTj{t')i--].HTAQ] 



(B.IO) 



S{t'i.-..,<_p,t=0) 
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We now recall that in the open system limit we take vp/L <^ l/|^o| ^ V with 
to — — oo. This implies that the e*** term in Ht^i essentially cuts off the lower 
limit of the time integrals at a time ^ much earlier than to- Thus, Eq. ( IB. 101 ) 
can be cast into an integral over {{A{tQ)B{t[, . . . ,t'^_p,t = 0)))o, where \t[ — 
to\,. . -M'n-p - ^oMO - tol ^ OO. Using Eq. we get 

{{A{to)B{t[, . . . , C,, t = 0)))o = {{A{to)))o{{Bit[, . . . , Cp, t = 0)))o, 

(B.ll) 



and hence 

{0)r 



tr [pn,/(O)0j 

tr[po] 



dtl 



to 



n—p 



dt' 



to 



n—p—1 



[[. . .[Oj{0),HTAtn^p)] ■ ■ ■ ,]HTAt[)]PO 

Inverting the nested commutator structure as previously, we get 



{0)n 



^ tr[pp,,(to)] 



^ "'to 



" dt' 

to 



p=i tr[po 

[ifT,/(Cp),[---,[^T,/(t'i),P0]]...]O/(0) 



(B.12) 



to 



dt\ 



(B.13) 



to 



(B.14) 



Finally, comparing the second term in the latter equation with Eq. (IB. II ) and noting 
that tr[pp^/(to)] = tr[pp /(O)], we arrive at the final recursion relation 



{0)n 



tr[Po] 



h tr[po 



(B.15) 



C Two derivations of Ud 



C.l Derivation ofud using tJie spectral representation 



Here, we provide additional details on the equivalence between the effective equi- 
librium approach and the Schwinger-Keldysh formalism and explicitly show the 
identity of the electron occupancy on the quantum dot as obtained in the Schwinger- 
Keldysh formalism and as obtained with the effective equilibrium approach. 
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We start by stating two useful identities involving Liouvillian superoperators. The 
first identity regards the imaginary-time Heisenberg representation of an operator 
O and is given by 



(C.l) 



The second identity is a small trick which transfers a Liouvillian superoperator 
from one side of an anticommutator to the other, 

{0„-1^04) = ({^0,.04). (C.2) 

where 2 is a c-number. 

We now set out to demonstrate the equivalence of the Schwinger-Keldysh expres- 
sion for the dot occupancy, 



/oo 
deAdM 
-00 



2 ^ j-00 
1 



2/ " V 2 

and the expression obtained in the effective equilibrium approach. 



(C.3) 



(C.4) 



It is crucial to recall the definition of the Green's functions in the effective equi- 
librium approach (see Section [3^ . The Fourier representation of the (imaginary) 
time ordered Green's function is given by 



iUn — C + Ly 



-a 



(C.5) 



and the retarded/advanced real-time Green's function in frequency space is ob- 
tained by 



Gret/adv/, ,\ 



1 



CO ^ C ±ir] 



(C.6) 



To prove the equivalence between Eqs. (IC.3I) and (|C.4I) . we start with the Keldysh 
expression and evaluate the spectral function 



Ad(e) = —lmG 



TT 



1 



^adv 



Gret 

dad'a 



{C.l) 
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With Eq. (IC.6I) we evaluate the retarded and advanced Green's functions involved, 



where, for simplicity, we are considering symmetric tunnel couplings and have set 
ti = t„i = t. Similarly, one finds for the advanced Green's function 



= ({^- + = T^({4,^«..}). (C.9) 

Utilizing these relations in the evaluation of the occupancy, we find 



^'^ = J-T7^ E (({4, i^ak.}) - h.c.) / (e, - a^) . (C.IO) 



47rztzy(0)Vl]^ 
This is equivalent to: 



111 gj'<JnO+ 



47rUzy(0)v^/3 lUJn - efc + 



a 



2 



47ritz/(0)v^/3 V \ I + ~ -^^^ "'^1/ \ [ «w„ - £ + £ 



47rztz/(0)v^/3 .ti^ V \ r ' + - A 

+ /Lt e-""°^ t 1 ^ n 

\\ zw„ + £ - £y Vfiefc + £ - z?7 '^j/ 

-Cafc^, f - ( ^ — ir—r 7^- 7r-—dl, 



iun — C + Cy'^"'^"''"'' \ I \\iujn — C + Cy y/U^k — + irj \ j j 

(C.ll) 



The first and third terms can be combined, and recalling Eq. (1191 ) they can be shown 
to cancel. 



11^/^ . ^ .\ 1 



(C.12) 
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since in the steady state we must satisfy Ji = Therefore, this results in: 



rid 



1 



1 



^o+ 



1 



.0+ 



iujn + C — Cy ek + C — if] 
1 



S d 

ZOOn — L + Ly ek — L + IT] 



47rz P ^-oo 



dei. 



e 



-dl 



1 



[ iuJn — C + Ly tk + L — ITj 



-4, 



iun - C + ir] " ek + C - C 



-dn 



Y 



(C.13) 



In the next step, we will make use of the identity 



oo 1 

de O = ±nvO, 

-oo e ± ir] 



(C.14) 



and we briefly outline its proof. It is convenient to use the spectral representation 
for this purpose and thus enumerate the set of eigenstates of if by (In)}. With this, 
one obtains 



de {m\ 0\n) = / de ' 

-oo e± iri J-oo e ± (e^ — Cn) =F 



m\0\n)de. (C.15) 



After shifting the integration variable e — )■ e ± (e^ — e„), we can further simplify 
the above expression to yield 



30 1 / 

de {m\0\n) = P.V. 

CO e =F zr] \ 



~ de 

-oo 6 



± in {m\0\n). 



(C.16) 



Noting that the principal value integral vanishes concludes the proof of Eq. (IC.14I ). 
With this identity in place, Eq. (IC.13I) can be brought into its final form 



1 ^ 1 

aijJnO- 1^ 



^0+ 



_ \ lUJn + C- Ly ] I 



-■ r I r 4^g^a 

lUJn — L + Ly 



(C.17) 



This demonstrates that the Schwinger-Keldysh approach and the effective equilib- 
rium formulation reproduce the same expression for the occupancy on the quantum 
dot. 
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C.2 Derivation ofud using the Keldysh approach 



The basic equation which will constitute the starting point of our proof is given by 
Eq. (15) in Ref. lfT4ll . The current in lead a = ±1 for spin projection a is given by 

X [G<At,h) + fa{e)GUt,ti)]'^, (C.18) 

where we have restricted ourselves to a single level on the dot having a spin in- 
dex a. The occupation number(r2fj^) of the dot in the state a obeys the following 
differential equation: 

TL = l^alaa- (C.19) 

at -e „ 

In steady state (i.e., t=0) we get 

dridAt) 



= 0. (C.20) 



dt 

This implies 

E £^^^1 / ^Im|e-"*T"(0,ti) [G<.(0,ti) + Ue)GUO,h)] | = 0. 

(C.21) 

Since in steady state there exists time translational invariance, we have ^^^(0, ti) = 
G^^{0-ti) and 6*^^(0,^1) = G'^^(O-ti). Further we can set 77 -> and to -00, 
in the sense defined by the open system limit, right from the beginning. This in turn 
implies that r"(0, ti) = = 7ri^(0)t^ = T/2. Let us focus on the first term in the 
expression above 



Using the fact [G<,(-ti)]* = -G<M 

(C.22) 
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Similarly for the second term we get 



r 

2 



2 TT 



/■=ff(£,4.) 



Combining the two terms and putting in Eq. (IC.21I) we get 



(C.23) 



(C.24) 



which implies that 



(C.25) 



This is identical to the expression obtained via the spectral representation in Ap- 
pendix [CT] 

It is important to mention at this point that the expression for the local charge in 
terms of the spectral function in Eq. (I30l) depends crucially on the fact that the 
system is in steady state. Unlike its equilibrium counterpart, the equation (T< (lo) = 
i[ar{ui) — (Ta(u;)]/''^'^'^) does not hold in general 117 8 1791 . Here (T<(u;), ar{uj) and 
aa{uj) denote the self energy contribution due to interactions to G^{uj), G^{uj) and 
G"-{u) respectively. Using Dyson's equations in Appendix IC. 21 it is simple to show 
that the weaker relation 



de 



a<(e) - i[ar{e) - cra{e)]f 



eff 



0, 



(C.26) 



[(e - e, - Re[a,(e)])2 + (r - Im[a,(e)])2j 
holds in steady state and leads to the equilbrium-like form of Ud given by Eq. (|30l) . 



D Second-order perturbation theory details 



In this appendix we will outline the steps involved in the evaluation of the retarded 
electron Green's function to second-order in Hi^t of Section 5.3. First we evaluate 
the relevant D^^^(c<;)'s involved in this computation. In the following discussion 
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we will use the abbreviated notation which was introduced after Eq. (1681) . i.e., I = 

{aiki),l = (aikf), gi = ^(e^J and gj = g{ekj). 



Vfi 1 - ei + z?7 



follows directly from the definition. Using Eq. (I63l) we get 



D^\^) = —,—-C,D^^\ (D.2) 

OJ + + IT] 

To compute this quantity we first evaluate the action of the superoperator Ci on ipj^. 

A^i. = -Ug^ E ^i52^3V^S°JI#l4°J - I E 9i4T^ ■ (D.3) 
Using this expression in Eq. (ID.2|) we obtain 

^ V ' 

i2 ^ w + ei - 62 - 63 + z?7 y 

Vfi \ i Z ^ W + 6i - 62 - €3 + ZT] 

-^E (D.4) 

2 ^ u; - 6i + ir] / 

The anticommutator of interest is then readily computed to give 

1 . . 



{D^^\u),dl} = Ug{u) 



(D.5) 



Similarly D^^\uj) follows from the use of Eq. (l63l) in Eq. (ID.4I) above 

(c^) = — (D.6) 
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We begin by computing 

U 



" 12ct \"/ 1234CT 



u 
2" 



12f7 



+ 2 



^ ^ V 1234CT 



(D.7) 



It is convenient to divide the action of the Liouvillian above into two parts by defin- 
ing 



, ' 12a- 

I)' E ^r^2*^3^4^£^f-^l#WlJ,^«' 

1234CT 



(D.8) 



We introduce the notation 



^i(!4)(^) 



1 



oj + C + ir] 

to separate the contributions to D^^^ (cu) which follow from the parts A and B. 



(D.9) 



The evaluation of the Green's function is likewise separated into two parts. First, 
we compute the contribution of the part A 



- / dei \gi 

71 J 



r/2 



1 



u — ei + 2ir uj + ei + 2ir 



fi 



eff 



. (D.IO) 



Similarly, we compute the the contribution of part B to the retarded electron Green's 
function of the dot G^^^^V (uj) 

dad' ^ ' 



3 (^) / A, / A2 ISil" |g2l'9(i^ + ti - t2)/r/t 



eff /eff 
2 



1 1 r r 

o , o r / ^1 



+ ei + 2ir 



eff 



. (D.ll) 
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Combining Eqs. (ID. 101) and (ID.l II) we finally obtain the second order correction to 
the retarded electron Green's function of the dot 



-yret(2) 

•A, 

t2 .m2 



-3 



r 



71 



oj + 3iT 

dei / de2 \gif {92^ 9 {^^ + ei - e2)fff2 
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E Deviations from the particle-hole symmetric point: Summary of results to 



We study the effect a deviation from the particle-hole symmetric point for the An- 
derson impurity 



a.k(7 



(E.l) 



Here 5^ measures the deviation from the particle-hole symmetric point. For a local 
magnetic field H on the dot, we have 



5„ = -a- 



H 



(E.2) 



The retarded electron Green's function of the dot now bears an explicit a depen- 
dence, i.e. 



9dA^) 



(E.3) 



We define the auxiliary function 



1 



9dA^} 



(E.4) 



Furthermore, the deviation of the occupancy of the dot level given by the spin pro- 
jection cr is denoted by 



(E.5) 
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The first-order contribution to the complete self-energy S^''^^) (defined in Eq. (|83l )) 
is frequency independent and is given by 

EW = [/5nd,. (E.6) 



In The second-order contribution to the complete self-energy consists of 3 parts. 
The contribution following from A in Appendix ID] gets modified to 



v(2) 



1 



u — El — 5^ + 6^^ + 2iT 



u + ei - 6^ - + 2iT j 
whereas the contribution which follows from B gets altered to 



(E.7) 



MB) 



r 



TT 



A I A \ |2 I |2 AeS) /•(eff) 

dei / dea \gia\ \g2-a\ Ji J2 



1 



1 



OJ + £2 -- Sa - 5-a + 2iT 
2 



Lo - €2 - 5„ + + 2ir 
+ 9d-A^ + €2 - ei) - Qd^A^ - €2 - ei)^ 

2 ^(eff) ^(eff) 

ydaV^ -r ci — C2 

1 



^) I J \92-aY fT'''f2'''9dA^ + ei - £2) - 9dA^)^'^d 



x[- [de, \g,^^f ff) -- [de, \g^.^\' ffl— 



(E.8) 



Finally we now get a second-order contribution from the term (^^D''J^\uj) , d^j 



(1) 



^(2) 



9dA^) 



TT 



de^n 



(eff) 



\9la\ 9l-a 



j_ \9l~cy\ 9la 

uj + ei — 5„ — + 2ir uj — ei + 2iT 



Sud^- 
(E.9) 



- (2) 

The contribution S^j.^-), corresponds to two concatenated Hartree-Fock diagrams, 
and is thus not contribute to the one-particle irreducible (IPI) diagrams, which de- 
fine the self-energy S'*^' (see Eq. (|82l) ). 



In other words, for the case of a local magnetic field on the dot, the self-energy to 
second-order is modified to 
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Fig. 1. Schematic setup (upper panel) and energy diagram (lower panel) of a generic 
quantum impurity model out of equilibrium. The system is at a temperature T, and 
$ = ;Ui — = eV denotes the voltage bias between the source and drain leads. The 
(bare) energy of the dot level is given by e^. The energy broadening of this level, given by 
the width V, is due to the tunnel coupling between dot and leads. 



(^) (b) 




Fig. 2. Here, (a) and (b) denote Hartree diagrams contributing to the current to first-order in 
the interaction strength. The Fock diagrams (c) and (d) are irrelevant due to the exclusion 
principle. The absence of exchange symmetry is depicted by distinguishing a particular end 
of the interaction Une by a dotted circle. 
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Fig. 3. Spectral function (obtained using the Bom approximation) for different values of 
the bias voltage for U/ {-nV) = 2, where we set T = 1. The hmits <^/r — ;> and <l>/r oo 
agree with the results obtained via the Schwinger-Keldysh scheme 1166 1671 . 




Fig. 4. Behavior of the spectral function (obtained using the Born approximation) for dif- 
ferent values of U, in the extreme limit $/r — )• oo, where we set F = 1. The results ai^e 
identical to the results obtained using the Schwinger-Keldysh formalism ll66l . 
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(o/r 

Fig. 5. Vaiiation of the spectral function with an applied magnetic field (H), obtained us- 
ing the second-order self-energy. Here T = I, ^ = 0, U = Svr and H is given in the 
units defined by gfii,, where g is the g-factor and /Xf, denotes the Bohr magneton. The inset 
shows the behavior of the Abriksov-Suhl resonance as a function of the magnetic field, for 
H = 0.8, 1.0 and 1.6. 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




5 10 15 20 



Fig. 6. The current-voltage curves for the Anderson model for U/T = 0.0, 1.0, 27r and 
Ait, where F = 1. The inset shows the behavior of the curves for low bias. In the limit 
$ — > the slope of the curves tend to 1, which corresponds to the value of the conductance 
quantum. 
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